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Quick chemical equilibration times of hadrons (specifically, pp , KK , AA, and QQ pairs) within a 
hadron gas are explained dynamically using Hagedorn states, which drive particles into equilibrium 
close to the critical temperature. Within this scheme, we use master equations and derive various 
analytical estimates for the chemical equilibration times. We compare our model to recent lattice 
results and find that for both T c — 176 MeV and T c — 196 MeV, the hadrons can reach chemical 
equilibrium almost immediately, well before the chemical freeze-out temperatures found in thermal 
fits for a hadron gas without Hagedorn states. Furthermore the ratios p/ir, K/n , A/n, and Q/tt 
match experimental values well in our dynamical scenario. 



I. INTRODUCTION 

(Anti-)strangeness enhancement was first observed at 
CERN-SPS energies by comparing anti-hyperons, multi- 
strange baryons, and kaons to pp-data. It was consid- 
ered a signature for quark gluon plasma (QGP) because, 
using binary strangeness production and exchange reac- 
tions, chemical equilibrium could not be reached within 
a standard hadron gas phase, i.e., the chemical equilibra- 
tion time was on the order of r ~ 100 — 1000 ^ whereas 
the lifetime of a fireball in the hadronic stages is only 
t « 4 — pp. It was then proposed that there exists 
a strong hint for QGP at SPS because strange quarks 
can be produced more abundantly by gluon fusion, which 
would account for strangeness enhancement following 
hadronization and rescattering of strange quarks. Later, 
however, multi-mesonic reactions were used to explain 
secondary production of p and anti-hyperons [21 |3J . At 
SPS they give a chemical equilibration time Tp « 1 — 
using an annihilation cross section of a p y w a p p w 50mb 
and a baryon density of pb ~ po to 2po, which is typical 
for evolving strongly interacting matter at SPS before 
chemical freeze-out. Therefore, the time scale is short 
enough to account for chemical equilibration within a 
cooling hadronic fireball at SPS. 

A problem arises when the same multi-mesonic reac- 
tions were employed in the hadron gas phase at RHIC 
temperatures where experiments again show that the 
particle abundances reach chemical equilibration close to 
the phase transition g]. At RHIC at T = 170 MeV, 
where a w 30mb and p e g sa p e ^ ss 0.04fm~ 3 , the equilib- 
rium rate for (anti-)baryon production is t 10—. 

Moreover, r « 10— was also obtained in Ref. [5] using 
a fluctuation-dissipation theorem. From hadron cascades 
a significant deviation was found from the chemically 
saturated strange (anti-) baryons yields in the 5% most 
central Au-Au collisions [DJ. These discrepancies suggest 
that hadrons are "born" into equilibrium, i.e., the system 
is already in a chemically frozen out state at the end of 
the phase transition [7J [S] . In order to circumvent such 



long time scales it was suggested that near T c there ex- 
ists an extra large particle density over/populated with pi- 
ons and kaons, which drive the baryons/ anti-baryons into 
equilibrium [DJ. But it is not clear how this overpopula- 
tion should appear, and how the subsequent population 
of (anti-)baryons would follow. Moreover, the overpopu- 
lated (anti-)baryons do not later disappear [ID]. There- 
fore, it was conjectured that Hagedorn resonances (heavy 
resonances near T c with an exponential mass spectrum) 
could account for the extra (anti-)baryons [TP ] ITT ] 112]. 
Hadrons can develop according to 

nn <^ HS <-> n' it + XX (1) 

where XX can be substituted with pp , KK , A A, or 
f2f2. Eq. (JT|) provides an efficient method for producing 
of XX pairs because of the large decay widths of the 
Hagedorn states. In Eq. (JTJ, n is the number of pions for 
the decay nir «-> HS and n' is the number of pions that 
a Hagedorn state will decay into when a XX is present. 
Since Hagedorn resonances are highly unstable, the phase 
space for multi-particle decays drastically increases when 
the mass increases. Therefore, the resonances catalyze 
rapid equilibration of XX near T c and die out moderately 
below T c [TTJ. 

Unlike in pure glue SU{3) gauge theory where the 
Polyakov loop is the order parameter for the deconfinc- 
ment transition (which is weakly first-order), the rapid 
crossover seen on lattice calculations involving dynam- 
ical fermions indicates that there is not a well defined 
order parameter that can distinguish the confined phase 
from the deconfincd phase. Because of this it is natu- 
ral to look for a hadronic mechanism for quick chemical 
equilibration near the phase transition. One such pos- 
sibility could be the inclusion of Hagedorn states. Re- 
cently, Hagedorn states have been shown to contribute 
to the physical description of a hadron gas close to T c . 
The inclusion of Hagedorn states leads to a low rj/s in 
the hadron gas phase [13], which nears the string the- 
ory bound rj/s = l/(47r) [14]. Calculations of the trace 
anomaly including Hagedorn states also fits recent lat- 
tice results well and correctly describe the minimum of 
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the speed of sound squared, c 2 s , near the phase transition 
found on the lattice [13]. Estimates for the bulk viscos- 
ity including Hagedorn states in the hadron gas phase 
indicate that the bulk viscosity, £/s, increases near T c , 
which agrees with the general analysis done in 15j. Fur- 
thermore, it has been shown [16] that Hagedorn states 
provide a better fit within a thermal model to the hadron 
yield particle ratios. Additionally, Hagedorn states pro- 
vide a mechanism to relate T c and T c ^ em , which then 
leads to the suggestion that a lower critical temperature 
could possibly be preferred, according to the thermal fits 

Previously, in Ref. [TT] we presented analytical results, 
which we will derive in detail here. Moreover, we saw that 
both the baryons and kaons equilibrated quickly within 
an expanding fireball. The initial saturation of pions, 
Hagedorn states, baryons, and kaons played no significant 
role in the ratios such as K/n and (B + Bj /it. 

Here we consider the effects of various initial condi- 
tions on the chemical freeze-out temperature and we find 
that while they play a small role on the total particle 
number, they still reproduce fast chemical equilibration 
times. Additionally, we assume lattice values of the criti- 
cal temperatures (T c = 176 MeV [T7] and T c = 196 MeV 
[151 IT!?] ) and find that chemical equilibrium abundances 
are still reached close to the temperature given by ther- 
mal fits (T « 160 MeV). 

This paper is structured in the following manner. In 
Section [IT] we discuss the details of our statistical model 
that calculates the chemical equilibrium values of the 
Hagedorn states and other hadrons. Furthermore in 
this section, fits are shown to thermodynamical prop- 
erties calculated in lattice QCD, which are used to de- 
termine the mass spectrum of the Hagedorn states and 
the rate equations are discussed in detail. In Section III 
we are able to extract the chemical equilibration time of 
an XX pair when the pions and Hagedorn states are 
held constant. In Section IV we derive an analytical 
result of the rate equations when we consider only the 
decay HS «-> nir. We then discuss the case of an ex- 
panding fireball and the results for the various XX pair 
production in Section [V] The production of particles 
will also be considered in Section I VI I We summarized 
and discussed our results in Section |VII[ In Appendix [X] 
we present some analytical and numerical results for the 
various equilibration stages in the hadron and Hagedorn 
states gas mixture. 



II. MODEL 

Hagedorn resonances have an exponentially growing 
mass spectrum |20j . Their large masses open up the 
phase space for multi-particle decays. Recent analysis 
involving Hagedorn states is given in [21]. Moreover, 
thoughts on observing Hagedorn states in experiments 
are given in [35] and their usage as a thermostat in [23J . 
Hagedorn states can also explain the phase transition 



above the critical temperature and, depending on the 
intrinsic parameters, the order of the phase transition 
[24 . For the following discussion, the overall density of 
Hagedorn states in our extended Hagedorn gas model are 
straightforwardly described by, 
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where M = 2 GeV and m 2 = 0.5 GeV. We note that 
in this work we consider only mesonic Hagedorn states 
with no net strangeness. The exponential in Eq. ^ 
arises from Hagedorn's original idea that there is an ex- 
ponentially growing mass spectrum. Thus, as Tjj is ap- 
proached, Hagedorn states become increasingly more rel- 
evant and heavier resonances "appear" . The factor in 
front of the exponential has various forms [5TJH3J. While 
the choice in this factor can vary, it was found in |21j that 
the present form gives lower values of Th, which more 
closely match the predicted lattice critical temperature 

[Hj[THi[ig. 

Returning to Eq. ([2]), its parameters (A, M, and Th) 
are dependent on the critical temperature. We assume 
that Th — T c , and then we consider the two different 
different lattice results for T c : T c = 196 MeV [TSJ QJ], 
which uses an almost physical pion mass, and T c = 176 
MeV [17]. Furthermore, we need to take into account the 
repulsive interactions and, therefore, we use the following 
volume corrections (as also seen in [131 1231 12"6] V 
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which ensure that the our model is thcrmodynamically 
consistent . Note that B is a free parameter that is based 
upon the idea of the MIT bag constant. 

In order to find the maximum Hagedorn state mass 
M and the "degeneracy" A, we fit our model to the 
thermodynamic properties of the lattice. In the RBC- 
Biclcfeld collaboration the thermodynamical properties 
are derived from the quantity e — 3p, the so called inter- 
action measure, which is what we fit in order to obtain 
the parameters for the Hagedorn states. Thus, we ob- 
tain T H = 196 MeV, A = O.bGeV 3 / 2 , M = 12 GeV, and 
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B = (340GW) 4 . The fit for the trace anomaly 6/T 4 is 
shown in Fig. [T] We also show the fit for the entropy den- 
sity in Fig. [2] Both fits are within the error of lattice and 
mimic the behavior of the lattice results. As discussed in 
[T3"] . a hadron resonance gas model with Hagedorn states 
uniquely fits the lattice data whereas a hadron resonance 
gas without Hagedorn states completely misses the be- 
havior. 
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FIG. 1: Comparison of trace anomaly to lattice QCD results 
from EES [19] where T c = 196 MeV. 
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FIG. 2: Comparison of entropy density to lattice QCD results 
from [181 OH where T c = 196 MeV. 



BMW calculates the thermodynamical properties sep- 
arately and, therefore, we fit only the energy density as 
shown in Fig. [3] From that we obtain Tjj = 176 MeV, 
A = 0.lGeV 3 ^M = 12 GeV, and B = (300GeV) 4 . We 
also show a comparison to the entropy density in Fig. 
[4] Our results with the inclusion of Hagedorn states are 
able to match lattice data near the critical temperature 
but do not match as well at lower temperatures in Fig. [T] 
and Fig. [2j 

Our idea is that these very massive Hagedorn states ex- 
ist, as pictured in Fig. [5] and are so large that they decay 
almost immediately into multiple pions and XX pairs. 
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FIG. 3: Comparison of energy density to lattice QCD results 
from Q7] where T c = 176 MeV. 
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FIG. 4: Comparison of entropy density to lattice QCD results 
from IT] where T c = 176 MeV. 



While it can be argued that Hagedorn states are more 
likely to decay into a pair of particles: a lighter Hage- 
dorn state and another particle, these reactions are so 
quick that we can consider the end results, which would 
be multiple particles (mostly pions). That being said, it 
would be possible to put Hagedorn states into a trans- 
port approach such as UrQMD [27 using binary reactions 
with possible cross-sections as described in [28 . We leave 
this as a challenge for the future. 

Moreover, we need to consider the back reactions of 
multiple particles combining to form a Hagedorn state in 
order to preserve detailed balance. Rate equations pro- 
vide us with a perfect tool for this because there is a loss 
and gain term that describe both the forward and back 
reactions. Moreover, the state of chemical equilibrium is 
a fixed point of the rate equations. The rate equations 
for the Hagedorn resonances N i: pions N^, and the XX 
pair N X x, respectively, are given by 
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FIG. 5: Hagedorn states decay into multiple pions and a XX 
pair. 
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resonance are I\ ff and 



The decay widths for the i th 
T iXX , the branching ratio is Bi >n (see below), and the 
average number of pions that each resonance will decay 
into is (n.i). The equilibrium values N eq are both temper- 
ature and chemical potential dependent. However, here 
we set /ib — 0. Eq. Q can also be rewritten in terms of 
fugacities (Aj, A^, and X xx ), which are found by dividing 
each total number by its respective equilibrium value, for 
example, A^ = -^ft, (as seen for the baryon anti-baryon 

pairs in [TT]). Additionally, a discrete spectrum of Hage- 
dorn states is considered, which is separated into mass 
bins of 100 MeV. Each bin is described by its own rate 
equation. 

The branching ratios, B ini are the probability that 
the i th Hagedorn state will decay into n pions. Since we 
are dealing with probabilities, then J2 n B*,n — 1 must 
always hold. In order to include a distribution for our 
branching ratios we assume that they follow a Gaussian 
distribution for the reaction HS «-> mr 
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/2tt 



(5) 



which has its peak centered at (m) and the width of 
the distribution is a 2 . Assuming a statistical, micro- 
canonical branching for the decay of Hagedorn states, 
we can take a linear fit to the average number of pions in 



Fig. 1 in Ref. [TO] (multiplying tt + by three to include all 
pions) to find (n,) such that (m) = 0.9+1.2^- is the av- 
erage pion number that each Hagedorn state decays into. 
Within the microcanonical model a Hagedorn state is 
defined by its mass and corresponding volume where the 
volume is taken as V = rrii/s. The mean energy density 
of a Hagedorn state is e (taken as e = 0.5 j^)- Fur- 
ther discussions regarding this can be found in [TUl l2"9"] . 
The width of the distribution is of = (0.5^) 2 . Both 

£ N trip ' 

our choice in (rii) and erf roughly match the canonical 
description in [57] . 

Furthermore, we have the condition that each Hage- 
dorn resonance must decay into at least 2 pions. Because 
of the nature of a Gaussian distribution there is a non- 
zero probability that a Hagedorn state can decay into less 
than 2 pions. Therefore, we calculate the percentage of 
the distribution that falls below 2 pions and redistribute 
that over n > 2 so that ^2 n B i n = 1. This in turn 
leads to a new (n^) and <jf, which we find by calculating 
( n i) = ^2n n Bi,n and of = (nf) — (rii) 2 . Thus, we after 
normalize for the cutoff n > 2, we have (rij) ~ 3 — 34 and 
of w 0.8 - 510. 

For the average number of pions when a XX pair is 
present, we again refer to the micro-canonical model in 

2~|[10l|29]. We use (n n ) but then readjust it to the average 
pion number according to Fig. 2 in Ref. |10j for when a 
.baryon anti-baryon pair is present (there the distribution 
is for a resonance of mass m = 4 GeV). Thus. 

(4) 

M (0.3 + 0.4^2-7. (6) 

where rrii is in GeV. In this paper we do not consider a 
distribution but rather only the average number of pions 
when a XX pair is present. We assume that (ni lX ) = 
( n i, P ) = (ni.k) = (»Ka) = ( n i,n) for when a kaon anti- 
kaon pair, AA, or f2f2 pair is present. Ideally, (n^k), 
(iii t \), and (n^n) should be derived separately and will 
be done in a future paper using a canonical model [3T] . 

We used a linear fit for the total decay width similar 
to that used in Ref. ■ The total decay width 



0.15m, - 0.0584 



(7) 



(Ti and m, in terms of GeV), which ranges from r, = 
250 MeV to 1800 MeV, is a linear fit extrapolated from 
the data in Ref. [33] • However, in Eq. Q the total decay 
width is separated into two parts: one for the reactions 
HS <-> nn, r i x , and one for the reaction in Eq. ([T]), 
T i xx , whereby Ti = r,- i7r + Tj xx . Then relative decay 
width r,j xx is the average number of XX in the system 
(X) multiplied by the total decay width Ti. Essentially, 
a fraction of the decay of the i th Hagedorn state goes into 
XX (set by the number of XX the i th Hagedorn state on 
average decays into) and the remainder goes into pions. 

We find (p) by linearly fitting the proton in Fig. 2 in 
Ref. [TD] so that 
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0.058 m, -0.10 
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where rnj is in GeV and (p) w 0.01 — 0.6. Thus, T^pp 
is between 3 and 1000 MeV. Clearly, I\ )jr is then r,* i7r = 
— (p)r i 7r . Analogously for the kaons, the decay width 
is T i K g — (K)Ti where 



K+ = 0.075 m t + 0.047 



(9) 



where nij is in GeV, which is also taken from Fig. 2 in 
Ref. [ID]. We find that (K) = 0.2 to 0.95 [TQlES]. Thus, 
T t,KK is between 50 and 1700 MeV. 

For A we use a canonical model assuming that the 
baryon number B = 0, the strangeness 5 = 0, and the 
electrical charge Q = in order to calculate the average 
lambda number. The results of this are shown in Fig. [6] 
We find that our (A) is lower than that from the micro- 
canonical ensemble in [10 , which is also shown in Fig. 
IB] This corresponds to a decay width of r ; aa = 3 — 250 
MeV. 

Furthermore, the average number of f2's is also shown 
in Fig. [6] from our canonical model again assuming that 
the baryon number B = 0, the strangeness 5 = 0, and 
the electrical charge Q = 0. In Fig. [6] we multiple (f2) 
in order to better view the results. The resulting decay 
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The equilibrium values are found using a statistical 
model [34 , which includes 104 particles from the the 
PDG [33 (only light and strange particles). As in [3"4"j . 
we also consider the effects of feeding (the contributions 
of higher lying resonances such as the p or uj resonances 
on the number of "pions" in our system, i.e., N% q in- 
cludes "all" the pions from resonances from the PDG 
|33j). Feeding is also considered for the protons, kaons, 
and lambdas. Additionally, throughout this paper our 
initial conditions are the various fugacities at to (at the 
point of the phase transition into the hadron gas phase) 



a = X n (t ) , ^ = \i(t ) ,and0 = Xxx^o) , (10) 
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FIG. 6: Average number of A's and f2's. The fi's are calcu- 
lated within our canonical ensemble and the A's are calculated 
both in our canonical ensemble and a micro-canonical ensem- 
ble. 



which are chosen by holding the contribution to the total 
entropy from the Hagedorn states and pions constant i.e. 



SHad(T Q ,a)V(t ) + s HS (T ,{3i)V(t Q ) 
= s H ad+Hs(T )V(t ) = const . 



(11) 



and the corresponding initial condition configurations we 
choose later can be seen in Tab. |TT] SHad(To, a) is the en- 
tropy density at the initial temperature, i.e., the critical 
temperature multiplied by our choice in a. Because the 
hadron resonance is dominated by pions we can assume 
that a represents the initial fraction of pions in equi- 
librium. shs{Tq, Pi) represents the entropy contribution 
from the Hagedorn states at T c multiplied by the initial 
fraction of Hagedorn states in equilibrium. We hold a 
constantly and then find the appropriate Pi . The volume 
expansion, V(t) is discussed in detail following section 
entitled 'Expanding Fireball'. 



III. CHEMICAL EQUILIBRATION TIME 
ESTIMATE 



As a starting point of our analysis, we first estimate 
the chemical equilibration time of the XX by looking at 
the fugacity of the XX rate equation, i.e., Eq. Q can 
be rewritten in terms of A as shown for BB in Eq. (3) in 
[TTj . when both the pions and Hagedorn states are held 
constant. The XX rate equation then becomes 
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which we can integrate 
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and Ajfjf(O) = <j). Substituting 



in a = 1 and Pi — 1 when the pions and Hagedorn states 
are in chemical equilibrium, we rederive Eq. (7) in Ref. 
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which is shown in Fig. [7] From Eq. ( 15 1 we sec that 
the time scale has an indirect dependence on the decay 
width. Since the decay width has a linear dependence 
on the mass, the time scale decreases when more Hage- 
dorn states are included. However, Nf q also decreases 
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FIG. 7: Comparison of the chemical equilibrium times for 
protons, kaons, and lambdas when a = 1 and Pi = 1 where 
T H = 176 MeV (top) and T H = 196 MeV (top). The gray 
band is the range of chemical equilibrium times for the Hage- 
dorn states (see Tab. [l| . 



with increasing mass so above a certain point very many 
Hagedorn states need to be included in order to see an 
effect in the time scale. Furthermore, the chemical equi- 
librium values have a dependence on the temperature, 
which makes the time scale shortest for the highest tem- 
peratures. 

In Fig. [8] and Fig. [9] we hold the Hagedorn states 
and pions and let the XX pairs reach chemical equi- 
librium. That means that in Eq. Q we set = 
and Ni = N^ q in the N X x equation. Fig. ^ shows the 
results for pp, KK, and A A, respectively, for Tjj = 176 
MeV and Fig. [9] shows the same for T H = 196 MeV. In 
all cases the temperature is held constant while the rate 
equations are solved over time. At T = T c all XX reach 
chemical equilibrium almost immediately (on the order 
of t < 0.2/m/c). As T is decreased the chemical equilib- 
rium time obviously increases, which is clear from Fig. 

ra 



Even as the temperature is lowered we still see quick 
chemical equilibrium times. For the pp and AA pairs at 
T = 0.9T C the chemical equilibrium time is still about 
t < 1/m/c. The KK pairs do have a slower chemi- 
cal equilibrium time due to their larger chemical equilib- 
rium abundances, which is directly related to the chem- 
ical equilibration time through Eq. (Q. This again rep- 
resents the main idea, which is the importance of po- 
tential Hagedorn states in understanding fast chemical 
equilibration of hadrons close and below T c . The Hage- 
dorn states increase dramatically in number close to the 
critical temperature and, thus, by its subsequent decay 
and re-population they will quickly produce the various 
hadronic particles. 

The equilibration of XX pairs then shown in Fig. [8] 
and Fig. [9] where the analytical result in Eq. ( fl3| matches 
the numerical result exactly. From Fig. [8] and Fig. [9] it 
can be seen that all XX pairs equilibrate quickly close 
to the critical temperature r < l^p. Clearly, though, 
as the temperature decreases the chemical equilibration 
time lengthens. However, at Tjj — 196 MeV chemical 
equilibrium is still reached quickly, r < 1^. 



IV. ANALYTICAL RESULTS: PIONS AND 
HAGEDORN RESONANCES 

While the chemical equilibration time derived in the 
previous section is a good estimate, it can only be strictly 
applied when the pions and Hagedorn states are assumed 
to stay in chemical equilibrium at a constant temperature 
(Fig. [8] and Fig. [9]). Otherwise, non- linear effects that 
appear when the pions and Hagedorn states are allowed 
to equilibrate appear. 

To understand the dynamics in more detail, we con- 
sider the simplified case when the Hagedorn resonances 
decay only into pions HS «-> mr, which gives 
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N eq I N* 
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Assuming that the pions and the Hagedorn states de- 
scribed in Eq. ( 16 1 are then allowed to equilibrate near 
T c in a static system, we are able to derive analyti- 
cal solutions, the derivation of which is shown in detail 
in Appendix |Aj For the analytical solutions we divide 
the chemical equilibration into three stages, the chemical 
equilibration times of which are shown in Tab.[T] The first 
stage (described by r° in Tab.[l| of the evolution is domi- 
nated by the chemical equilibration of the pions when the 
pions are still far away from their chemical equilibrium 
values. After the pions are close to chemical equilibrium, 
new dynamics take over, which are described by r n in 
Tab. H] and Fig. [To] 

In both Stage 1 and 2 the equilibration of the Hage- 
dorn states is set by the dynamics of the pions. Finally, 
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FIG. 8: Graph of the number of proton anti-proton pairs, 
kaon anti-kaon pairs, and lambda anti-lambda pairs when 
both the resonances and pions are held in equilibrium for 
T H = 176 MeV. 

in Stage 3 the pions, which are already almost in chemi- 
cal equilibrium, reach a quasi-equilibrium state with the 
Hagedorn states. Quasi-equilibrium is reached when at 
least one species of Hagedorn states has succeeded it's 
chemical equilibrium time scale determined from the in- 
verse of its decay width, i.e., Tj = 1/IV Since the heav- 
iest Hagedorn states have the shortest Tj's, then quasi- 
equilibrium is reached when Tj of the heaviest Hagedorn 
state is surpassed. During this stage non-linear affects 
take over and, thus, a longer time scale, t^ e , is seen. 




t [fm/c] 

FIG. 9: Graph of the number of proton anti-proton pairs, 
kaon anti-kaon pairs, and lambda anti-lambda when both 
the resonances and pions are held in equilibrium for Th = 
196 MeV. 

While this time scale may appear long, both the pions 
and Hagedorn states are so close to chemical equilibrium 
that they are within roughly 10% (depending on the ini- 
tial conditions) or less of their chemical equilibrium val- 
ues before quasi-equilibrium is even reached. The de- 
tailed calculations are shown in the Appendix. 

Therefore, the most important chemical equilibration 
time is then that from the pions in Stage 1, i.e., The 
time scale from Stage 2 is so short that it is not of much 
importance. Additionally, by the time that Stage 3 is 
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M; 



0.8 J 



0.95^ 



2Gev Mi2Gev chemical equilibrium, while the lighter Hagedorn states 
are still only moderately close to their chemical equilib- 
rium values. Therefore, the Hagedorn states and pions 
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TABLE I: Chemical equilibration times from analytical esti- 
mates where QE is quasi-equilibrium at 95% of each respective 



are able to be roughly in chemical equilibrium on the or- 
der of < 1 — according to our analytical solution when 
held at a constant temperature. 

This also applies to the KK reaction HS <-» nir + KK 
as shown in the bottom panel of Fig. [31] in Appendix 
jAj The time scale for the pions and Hagedorn states are 
slightly longer when the KK pairs are present. The same 
goes for the estimated chemical equilibration time of the 
KK pairs in the previous section, t k ^. 
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V. EXPANDING FIREBALL 

In order to include the cooling of the fireball we need 
to find a relationship between the temperature and the 
time, i.e., T(t). To do this we apply a Bjorken expansion 
for which the total entropy is held constant 



const. = s(T)V(t) 



S„ f dN v 



dy 



dy. 



(17) 



where s(T) is the entropy density of the hadron gas with 
volume corrections. 

The total number of pions in the 5% most central col- 
lisions, can be found from experimental results in 
[35] . There they found the phase-space yields for the pi- 
ons 7r + (292.0) and 7r~ (290.9) using a Gaussian fit for 
yields as a function of the rapidity where we used 
the rapidity range y — ±0.5. We then assumed that the 
number of 7r°'s were also in that same range and took 
the average of the two to find 291.5. Thus, our total 
pion number is N„< = f°* 5 ^dy = 874. While 
for a gas of non-interacting Bose gas of massless pions 
Sn/Nn = 3.6, we do have a mass for a our pions, so we 
must adjust S^/N^ accordingly. In [33] it was shown 
that when the pions have a mass the ratio changes and, 
therefore, the entropy per pion is close to S^/N^ rs 5.5. 
The actual S^/N^ in our model is shown in Fig. 
Stt/N^ ss 6, which is only slightly higher. 

The effective volume at mid-rapidity can be 
parametrized as a function of time. We do this by 
using a Bjorken expansion and including accelerating 
radial flow. The volume term is then 
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where 



FIG. 10: Comparison of the chemical equilibration times of 
the pions to the total chemical equilibration time for Th = 
176 MeV (top) and T H = 196 MeV (bottom). 



reached both the pions and Hagedorn states are essen- 
tially in chemical equilibrium and, therefore, the non- 
linear affects do not play a large role in the overall chem- 
ical equilibration time. One can see this more clearly in 
the top panel of Fig. [31] in Appendix [A] where the pi- 
ons and heavier Hagedorn states are extremely close to 



F(/> = 7T <-/ ( r a + v (t - t Q ) + ^a {t - t ) 2 



(18) 



where the initial radius is ro(io) = 7.1 fm for Tjj = 196 
and the corresponding ig 196 ^ « 2/m/c. For Th = 176 
we allow for a longer expansion before the hadron gas 
phase is reached and, thus, calculate the appropriate 
from the expansion starting at Tjj = 196, which 



(176) 



is t, 



(176) 



4/m/c (there is a slightly variation dependent 
on the choice of Vq and et ). The T(t) relation is shown 



in Fig. 12 which has almost no effect on the results as 
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FIG. 11: Entropy per pion for a hadron gas in chemical equi- 
librium within the fireball ansatz. 
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FIG. 13: Ratio of the entropy of the Hagedorn states to the 
total entropy. 



seen later on in Fig. 16 and Fig. 17 Therefore, we choose 
vq = 0.5 and ao = 0.025 for the remainder of this paper. 
The relation depicted allows to translate the later shown 
figures labeled by the effective global temperature of the 
evolving system directly into the evolving system time. 
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FIG. 12: The temperature-time relationship is directly linked 



to the average transversel velocity chosen in Eq. ( 18 1 within 
the fireball model ansatz. 



contribution, we arrive at our effective number of pions 



AT 



,xx 



N~ 



[(1 - W) te) + (XiKni,*)] (19) 



which arc shown in Fig. [14] In Fig. [14] we see that after 
the inclusion of the effective pion numbers that the num- 
ber of pions only decreases with decreasing temperature. 
Furthermore, in Fig. [14] the total number of Hagedorn 
states, ^ N^ q is also shown. While there are by far 
fewer Hagedorn states present than pions, we see that 
they are important because of their large contribution to 
the entropy density as shown in Fig. [13] The reason that 
the effective number of pions increase close to T c is due 
to the large number of pions that the heavy Hagedorn 
states decay into. If (m) was smaller or no longer lin- 
ear than it could be possible that the effective number of 
pions would remain constant. 

Moreover, it is useful to consider the effective number 
of XX pairs 



N 



xx 



N 



xx 



E 



Ni(Xi 



(20) 



Because the volume expansion depends on the entropy 



according to Eq. (17) and the Hagedorn resonances con- 



tribute strongly to the entropy only close to the critical 
temperature (see Fig. 13 1, the equilibrium values actu- 



ally decrease with increasing temperature close to T c for 
the hadrons as seen in Fig. [14] and Fig. [15] One can 
clearly see from Fig. [13] that the Hagedorn states con- 
tribute strongly close to T c down to about 80% of T c . 

Therefore, one has to include the potential contribu- 
tion of the Hagedorn resonances to the pions as in the 
case of standard hadronic resonances, e.g. a p-meson de- 
cays dominantly into two pions and, thus, accounts for 
them by a factor two. This is similar to what was done in 
Appendix [A] in Eq. ( A7 ) . Including the Hagedorn state 



because Hagedorn states also contribute strongly to the 
XX pairs close to T c as seen in Fig. 15 Again we see 



that only the effective number of XX pairs have consis- 
tent decreasing behaviour with decreasing temperature 
whereas without the Hagedorn state contributions we see 
a decrease close to T c . 

Along with the expansion we also must solve these rate 
equations, Eq. Q, numerically (3S]. We start with var- 
ious initial conditions, as mentioned previously, that are 
described by a, Pi, and </> (see table II). The initial tem- 
perature is the respective critical temperature and we end 
the expansion at T — 110 MeV, a global kinetic freezeout 
temperature. 

For the remainder of this paper we include only results 
for an expanding fireball, which are solved numerically. 
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FIG. 14: Comparison of the effective pion numbers when 
T„ = 176 MeV (top) or T„ = 196 MeV (bottom). 
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TABLE II: Initial condition configurations, recalling Eq. (FLO 



As an initial test we hold both the pions and Hagedorn 
states in chemical equilibrium and allow just XX to equi- 
librate as seen in Fig. [16] and Fig. [17] The black solid line 
in each graph is the chemical equilibrium abundances 
and the colored lines are the dynamical calculations for 
various expansions that follow the T(t) shown in Fig. 
T2] We see that regardless of our volume expansion they 
all quickly approach equilibrium. In Fig. 16 and Fig. 
17 the XX all reach chemical equilibrium almost im- 



mediately, well before 0.9T C the chemnical equilibration 
time is < 1 — . The only exception is the KK pairs for 
Th = 176 MeV. However, we see later on that the K/ir 
ratio matches the data. 

More interestingly, we consider the case when the pi- 
ons, Hagedorn states, and XX all arc allowed to chemical 
equilibrate. We then vary the initial conditions and ob- 
serve their effects. The results for pp pairs are shown 
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FIG. 15: Comparison of the total number of XX and their 
effective numbers when Th = 176 MeV (top) or Th = 196 
MeV (bottom). 



in Fig. [18] and Fig. [19] In Fig. [18] and Fig. [19] we show 
the evolution of both the pp pairs and the pions for the 
reaction nir <-> HS <-> nir + XX. Note that in all the 
following figures the effective numbers are shown so that 
the contribution of the Hagedorn states is included. 

One can see that the chemical equilibration time does 
depend slightly on our choice of j3i, i.e., a larger (3i means 
a quicker chemical equilibration time. For instance, if 
the Hagedorn states were overpopulated coming out of 
the QGP phase than chemical equilibrium times would 
be slightly shorter. However, even when the Hagedorn 
resonances start underpopulated the pp pairs are able 
to reach chemical equilibrium immediately Addition- 
ally, when the pp pairs start at about half their chemical 
equilibrium values, it only helps the pp pairs to reach 
equilibrium at a slightly higher temperature (on the or- 
der of a couple of MeV). Additionally, we see a greater 
dependence on for Th = 176 MeV than for Th = 176 
MeV. Throughout the evolution we see from the pions 
that they remain roughly in chemical equilibrium. Thus, 
our initial analytical approximation appears reasonable. 

In Fig. [20] the ratio of protons's to 7r's is shown. We 
also compare our results to that of experimental data. 
We see that for Th = 176 MeV that our results enter 
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FIG. 16: Results for the pp, KK, and AA when the pions and FI G. 17: Results for the pp, KK, and AA when the pions and 

Hagedorn resonances are held in equilibrium for T H = 176 Hagedorn resonances are held in equilibrium for T H = 196 

MeV. MeV - 



the band of experimental data before T — 170 MeV and 
remain there throughout the entire expansion regardless 
of the initial conditions. However, for Tjj = 176 MeV 
the results are slightly different. In this case, the ratios 
match the experimental data early on at around T — 
190 MeV. However, they become briefly overpopulated 
around T = 160 — 170 MeV but then quickly return to the 
experimental values, except for the case when we have the 
initial conditions such that the pions are overpopulated. 
This could imply that there are a few too many Hagedorn 



states and a fit for the Hagedorn states with a lower A 
(degeneracy of the Hagedorn states) may produce better 
results. 

As with the protons, the total number of kaons are 
also slightly dependent on our chosen initial conditions, 
more specifically, our choice in /3j. In Fig. [21] and Fig. [22] 
the temperature of the evolving system after the phase 
transition at which chemical equilibrium among standard 
hadrons is basically reached and maintained is between 
T = 160 — 170 for T H = 176 MeV and they have also 
already reached chemical equilibrium by T — 170 for 
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FIG. 18: Results for the p's and pions with various initial 
conditions for T H = 176 MeV. 
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FIG. 19: Results for the p's and pions with various initial 
conditions for T H = 196 MeV.. 



Th = 196 MeV, below which the Hagedorn states ba- 
sically die out. The one exception is when the Hagedorn 
states begin underpopulated i.e. that /J; < 1. In this 
case, the kaon pairs take longer to reach chemical equi- 
librium. However, when we look at K/t: in Fig. 23 lower 
Pi actually fits the data better. 

Moreover, the pions again remain roughly at chemical 
equilibrium throughout the expansion as seen in Fig. pT] 



26 In both cases 



and Fig. 122 
in Figs. fl& 



While the pion graphs look roughly similar 
22~1 they are not. The difference is how the 



pions are affected in the presence of a pp pair compared 
to a decay that includes a kaon anti-kaon pair. 

In Fig. [23] the ratio of kaons to pions is shown for 
T H = 176 MeV and for T H = 196 MeV. For T H = 176 
MeV our results are roughly at the upper edge of the 
experimental values. However, for Th = 196 MeV our 
results are slightly higher than the experimental values. 
Although, the results at T — 110 MeV are almost exactly 
those of the uppermost experimental data point. 

We can also observe the affects of the expansion on the 
AA pairs as seen in Fig. 24 and Fig. 25 We see that both 
reach the experimental values almost immediately (T > 
170 for T H = 176 MeV and around T = 190 for T H = 196 
MeV) . The one exception is again for an underpopulation 
of Hagedorn states, which reaches chemical equilibrium 
at T « 165 for T H = 176 MeV and already by T = 170 



for T H = 196 MeV). 

The ratio of A/ V's is shown in Fig. 
the A/7r's match the experimental values extremely well. 
For Tjj — 176 MeV our results reach the equilibrium 
values at T ss 170 MeV and for T H = 196 MeV the ex- 
perimental values are reached already by T w 170 MeV. 

A summary graph of all our results is shown in Fig. 
[27] The gray error bars cover the range of error for the 
experimental data points from both STAR and PHENIX. 
The points show the range in values for the various initial 
conditions at T = 110 MeV. We see in our graph that 
our freezeout results match the experimental data well. 

What the graphs in Figs. [I8p6| show us is that a dy- 
namical scenario is able to explain chemical equilibration 
values that appear in thermal fits by T = 160 MeV. In 
general, Tjj = 176 MeV and Tjj = 196 give chemical 
freeze-out values in the range between T = 160 — 170 
MeV. These results agree well with the chemical freeze- 
out temperature found in |16j . 

Moreover, the initial conditions have little effect on 
the ratios and give a range in the chemical equilibrium 
temperature of about ~ 5 MeV, which implies that infor- 
mation from the QGP regarding multiplicities is washed 
out due to the rapid dynamics of Hagedorn states. Lower 
Pi does slow the chemical equilibrium time slightly. How- 
ever, as seen in Fig. [27] they still fit well within the ex- 
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FIG. 20: Results for the ratio of p's with various initial condi- 
tions. Note that for STAR p7r~O.ll and pir' = 0.082. Along 
the top axis of each graph the corresponding time is shown in 
fm/c. 



The first decay channel of a mesonic non-strange Hage- 
dorn state we can implement straightforwardly with our 
model by employing the canonical branching ratio via 
Fig. [6] The results are shown in Fig. [28] for T H = 176 
MeV, in Fig.|29|for T H = 196 MeV, and the Sl/V ratio is 
shown in Fig. 



30 



We are able to find the average num- 
ber of Si's from [Bl, as seen in Fig. 6j We see that, using 
only the first reaction, we are still impressively able to 
adequately populate the S1S1 pairs so that they roughly 
match the experimental data. On the other hand, from 
Fig. [7] we see that for the SI particle the equilibration 
time are short only very close to T c . The scenario is 
thus more delicate. If one would take eg one half, or one 
fourth, respectively, of the decay width of that of eq. [7] 
the total production of SI is not sufficient up to 25 %, or 
up to 50%, respectively, to meet the experimental yield 
(the other ratios are not significantly affected by such a 
change of the decay width) . 

In a future work, it would be interesting to observe the 
other decay channels as given in Eq. [2T]and advertised in 
[TP] . The second reaction includes a mesonic, three times 
strange Hagedorn state whereas the third decay channel 
includes a baryonic, strange Hagedorn state. Both states 
are much more likely to directly decay into a SI. These 
are, admittedly, exotic states, but should also occur in the 
spirit of Hagedorn states. In order to observe these decay 
channels a method, e.g. a microscopic quark model, must 
be found to find the appropriate Hagedorn spectrum for 
strange mesonic/baryonic Hagedorn states. 



VII. CONCLUSIONS 



pcrimental values. Furthermore, in |llj we showed the 
the initial condition play pretty much no roll whatsoever 
in the ratios of K/tt + and (B + B)/tt + . Thus, strength- 
ening our argument that the dynamics are washed out 
following the QGP. 

While the variance in the chemical equilibration time 
arising from the initial conditions may seem contradic- 
tory to the K/ir + and (B + B)/tt + ratios in [TT] , it can 
be explained with the pion populations. In Figs. [T8]|22| 
quicker chemical equilibration times and, thus, larger to- 
tal baryon/kaon numbers translated into a larger number 
of pions in the system. Thus, the K/tt + and (B + B)/ir + 
ratios do not depend on the initial conditions. 

VI. PRODUCTION OF Slfi 

We can also use our model to investigate the possibility 
of Si's. In [TO], they discussed the possibility of Si's being 
produced from the following decay channels: 

HS <-> S1S1 + AT 
HS(sssqqq) <-> fl + B + X 



In this paper we found that hadronic matter, at RHIC 
or SPS energies, can reach chemical equilibrium within 
a dynamical scenario using Hagedorn states close to the 
critical temperature. These states were able to produce 
quick chemical equilibration times in (anti-)proton, (anti- 
)kaons, and (anti-)lambdas close to the critical temper- 
ature due to their strong increase in their abundancy. 
The existence of such a mixture of standard hadrons 
and Hagedorn states just below the phase transition can 
explain dynamically the chemical equilibration of the 
hadronic species at around temperatures of 160 MeV to 
170 MeV as seen within the thermal models. 

From our analytical results we found that the chemi- 
cal equilibration time depends on the temperature, decay 
widths, and branching ratios, but not the initial condi- 
tions. While this changes slightly when an expanding 
fireball is considered, the initial condition still only play 
a small role and only minimally affect the 'freeze-out' 
temperature at which chemical equilibrium is reached. 
This demonstrates that regardless of the population of 
hadrons coming out of the QGP phase, the initial con- 
ditions are washed out and everything can reach abun- 
dances which correspond to those of chemical equilibrium 
by the chemical freezeout temperatures found in [16 . 
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FIG. 21: Results for the K's and pions with various initial 
conditions for T H = 176 MeV. 
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FIG. 22: Results for the K's and pions with various initial 
conditions for T H = 196 MeV.. 



Moreover, from our previous paper [TT] we showed that 
particle ratios (K/tt + and (B + B)/ir + ) are not affected 
by the initial conditions and here we showed that p/tt, 
K/tt , K/tt and also il/7r match the experimental val- 
ues regardless of the initial conditions. Specially, Fig. [27] 
demonstrates this nicely and summarizes our findings: 
Regardless of the initial conditions, our dynamical sce- 
nario can match experimental data. We do find, however, 
that Th — 196 fits within the experimental data box for 
K/tt whereas Th = 176 is slightly above. This appears 
to reconfirm the findings in [IB] . 

Our results imply that both lattice temperature can 
ensure that the hadrons reach their chemical equilibrium 
values by T = 160 - 170 MeV. Although the ratios for 
Th = 176 do fit the data somewhat better, both math 
the experimental values reasonably well. This implies 
that independent of the critical temperature the hadrons 
are able to reach chemical freeze-out. 

We see sufficiently short time scales for the chemical 
equilibrium of hadrons. The protons, kaons, and lambdas 
reach chemical equilibrium on the order of At w 1 — 2^. 
Moreover, Hagedorn states states provide a very efficient 
way for incorporating multi-hadronic interactions (with 
parton rearrangements). 

In an upcoming paper we will use a canonical model 
to derive all the branching ratios included in our cal- 
culations. We can then look at reactions that include 



a mixture of strange and non-strange baryons (for in- 
stance, p + A) and multi-strange baryons. However, con- 
sidering that our initial results produce quick chemical 
equilibration times for the baryons, kaons, and lambdas, 
it is reasonable to believe that this will occur for mixed 
reactions and multi-strange baryons as well. In addition, 
the machinery of standard hadronic reactions, i.e. binary 
scattering processes and resonance production processes, 
help also to equilibrate the various hadronic degrees of 
freedom. Still, our work indicates that the population 
and repopulation of potential Hagedorn states close to 
phase boundary can be the key source for a dynamical 
understanding of generating and chemically equilibrating 
the standard and measured hadrons. 
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APPENDIX A: APPENDIX: ANALYTICAL 
SOLUTIONS OF VARIOUS EQUILIBARTION 
PROCESSES 

If our initial conditions are such that both the pions 
and Hagedorn states begin far out of chemical equilib- 
rium, we can find an analytical solution by subdividing 
the analysis into three distinct stages. Initially, during 
stage 1 the pions are underpopulated such that we can 
say that they approximately begin at a f=a (we can 
also start the pions above zero and the approximation 
works well). Because the pions reach chemical equilib- 
rium much quicker than the Hagedorn states due to all 
the Hagedorn states decaying quickly into pions, then we 
can make the approximation that the Hagedorn states 
are held at their initial value of One can see this 
from the difference in the time scales from Tab. U where 
t s ; > t" and Tj > t„. Since a w we let A£ ~ 0, then 
substituting this into Eq. (16 I we obtain 



N eq 

^ t -Ti jy-eq Pi ( n i) ! 



(Al) 



which is the fugacity of the pions in stage 1 and gives 



Again using the approximation a 



equation in Eq. (16), with the solution 



and substituting Eq. ( Al I into the Hagedorn state rate 



A,; = 



1 - 



-t 

Ti 



(A2) 



Substituting x — — £ into the integral in Eq. (A2|, 
expanding the exponential inside the integral so e h = 
E^Lo and integrating over £, provides us with the 
fugacity of the Hagedorn states in stage 1 



<n*> 



1 -e 



3=0 



(A3) 



Therefore, Eq. ( |A1[ ) and Eq. (A3 1 describe the behaviour 
of the pions and Hagedorn states during the initial stage 
of the evolution towards chemical equilibrium. They are 
then compared to the numerical results in Fig. |31| 

As the pions near equilibrium our approximation of 
Att rj no longer holds and we switch to stage 2 where 
we assume A x ~ 1 at time t\. Here t\ is a time when 
the pions are almost in chemical equilibrium, which is 
normally taken when the pions reach about ~ 95% of 
their chemical equilibrium value. Returning to the pion 
equation in Eq. ( 16 1, we can substitute in = 1 — e and 
use the approximation (1 — e)" « 1 — ne 



N eq 

= -E r *}vH (A - 1)<r 



(A4) 



Additionally, we substituted in (3i for \ as an approxima- 
tion since the Hagedorn states do not change significantly 
in Stage 1 (the majority of the evolution is done by the 
pions). Recall that j3i = Xi(t = 0) and it is a constant. 
In its present form, Eq. (A4| can be integrated. We also 



define e(ii) = 1 — rj where r\ is close to 1 (77 is the mea- 
surement of how close the pions are to their equilibrium 
value when we switch from Stage 1 to Stage 2). Then, 
after integration 



At 



I + 7 — (I + 7- 77) e 



(A5) 
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FIG. 24: Results for the A's and pions with various initial 
conditions for T H = 176 MeV. 
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FIG. 25: Results for the A's and pions with various initial 
conditions for T H = 196 MeV.. 



where 7 



E<r 4 JV 4 e «(A-i)(rn) 



and 



N' q 



Analogously to stage 1, we substitute the pion equation, 
i.e., Eq. (A5l into the Hagedorn resonance equation in 



Eq. (16) and integrate 



de 



(m) (1 + 7 - r)) e 



(A6) 



= Ui-l + (n,)7 + (^^t) ( n *> (l + l-v) and 
LOi. Thus, our equations for the evolution of the 



where d 
Ai(*i) = 

pions and Hagedorn states are Eq. ( A5 1 and Eq. ( A6 1 , 
respectively. As with stage 1, the evolution equation for 
the Hagedorn states is dictated by that of the pions. 

Stage 3 i.e. quasi-equilibrium begins once the pions 
and at least one species of Hagedorn resonances (t tg&v 
is the shortest chemical equilibration time) has surpassed 
its equilibration time (t„ and tj, respectively). To under- 
stand quasi-equilibrium we must use the effective pion 
number 



N w = N„ + 



(A7) 



of both the pions and resonances. The effective pion 
number essentially includes the number of effective pions 
that each Hagedorn state could decay into. Thus, we 
start by taking the derivative of Eq. ( A7 1 in terms of its 
fugacity 



1 



eq 



N- 



eq 



(A8) 



Once again we make the substitution = 1 — e so that 



1 V 7T 



(A9) 



where erf = — (rii) 2 in the Gaussian distribution of 
our branching ratios. To relate e and e we return to Eq. 
(A7l and separate A^ into a sum over the resonances in 



because we need a variable that can observe the effects quasi-equilibrium and one over the "freely" equilibrating 
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FIG. 26: Results for the ratio of A/n's with various initial 
conditions. Note that for STAR A/tt~ = 0.54 and A/tt~ = 
0.41. Along the top axis of each graph the corresponding time 
is shown in fm/c. 
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FIG. 27: Plot of the various ratios including all initial condi- 
tions defined in Tab. [Ill The points show the ratios at T = 110 
MeV for the various initial conditions (circles are for Th = 176 
MeV and diamonds are for Th = 196 MeV). The experimental 
results for STAR and PHENIX are shown by the gray error 
bars. 



resonances 



1 



N^ + Y d Nt q {n i )\ i + 

QE 



free 



(A10) 

Since the pions reach quasi-equilibrium first, i.e., 7v < Tj 



near T c , we set the 7r rate equation in Eq. ( 16 1 equal to 
zero, which gives A 4 w -^-y J2 n B i,nnX%, so 



i 



Eq. (All ) then has the form A, « 1 - e where 



(All) 



(A12) 

We can then solve for e in Eq. i1A12|) and substitute e into 



Eq. (A9l, which in turn can be integrated. This leads us 
to the solution 



■ Y N t q (ni)e~^ [ e^Xi{x)dx. (A13) 

free 

where j stands for the latest resonance to reach chemical 
equilibrium at that point in time and 



Y,QE N i 9 ( n i) 



(A14) 



is the quasi-equilibrium time. Clearly, once all the Hage- 
dorn states have reached chemical equilibrium than j 
symbolizes the resonance of M = 2 GeV, since it is the 
slowest Hagedorn state to equilibrate. The sums over 
"free" is the sum over the Hagedorn states that have not 
yet surpassed their respective chemical equilibrium time, 
Tj. Once T2GeV is reached those sums equal zero. There- 
fore, after T 2 GeV au that remains is 



C2GeV e 



: '2< ; \ ' 

OE 



(A15) 



where t^ e is shown in Tab. [T] Finally, we rewrite Eq. 
(A15l in terms of the pion evolution equation 



A T = 1 — (1 — K)e 



'2a, r 



(A16) 



where k = X n {T2Gev)- 

Because the resonance equation depends on the pop- 
ulation of the pions we substitute Eq. ( A16 1 into the 



Hagedorn resonance rate equation in Eq. (16), assuming 
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FIG. 28: Results for the Si's and pions with various initial 
conditions for T H = 176 MeV. 




120 130 140 150 160 170 ISO 
T [MeV] 




140 160 ISO 
T [MeV] 



200 



FIG. 29: Results for the fi's and pions with various initial 
conditions for T H = 196 MeV.. 



the pions are near equilibrium (i.e., we use the approxi- 
mation A = 1 — e and (1 — e) n « 1 — ne) 



A; 



1 - 



1 



T QE -T 

'TV I 'I 



(rii)(l - k) e 



t QE - T 

I 7T I 1 



(ni)(l - n)e 



(A17) 



where &i = \i{r2Gev)- Thus, for stage 3 the population 
equations for the pions and the Hagedorn states are Eq. 
(A16l and Eq. (A17l so long as t > T 2 GeV- 



Fig. 31 reveals a remarkable close fit with our numer- 
175 MeV i.e. T < T H . Thus, the 



ical results for T 



quasi-chemical equilibrium time, t® e , depends only on 
Ti, (rii}, af, and N eq , which is temperature dependent, 
but not on our initial conditions. As mentioned in the 
text, though, t® e includes many non-linear affects that 
only occur close to the chemical equilibrium. Thus, the 
more appropriate time scale is t® in order to describe the 
dynamics. 



We also see from Fig. 31 that when KK pairs are in- 
cluded that the pions and Hagedorn resonances equili- 
brate in roughly the same amount of time, which implies 
that our analytical solution can still be approximately 
applied when KK pairs are present. 
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